Low Energy Quasiparticle Excitation in the Vortex State of Borocarbide 

Superconductor YNi 2 B 2 C 
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We have measured the low temperature heat capacity C p and microwave surface impedance Z s in 
the vortex state of YM2B2C. In contrast to conventional s-wave superconductors, C p shows a nearly 
x/TT-dependence. This %/77-dependence persists even after the introduction of the columnar defects 
which change the electronic structure of the vortex core regime dramatically and strongly disturb the 
regular vortex lattice. On the other hand, flux flow resistivity obtained from Z s is nearly proportional 
to H . Taken together, these results indicate that the vortex state of YNi2B2C is fundamentally 
different from the conventional s-wave counterparts, in that the delocalized quasiparticle states 
around the vortex core play a much more important role, similar to d-wave superconductors. 
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The borocarbide superconductors LixN^E^C where 
Ln=(Y, Lu, Tm, Er, Ho and Dy) exhibit a rich variety of 
interesting physical properties. In particular, the occur- 
rence of superconductivity at elevated temperatures |Q, 
the competition and coexistence of the antifcrromagnetic 
ordering and superconductivity g], and the transition 
between a triangular and square vortex lattice [g[ have 
attracted much attention. In spite of extensive studies 
on these subjects, however, one of the most fundamen- 
tal properties of the superconducting state, namely the 
quasiparticle (QP) structure in the vortex state, is still 
controversial. In fact, recent measurements of low tem- 
perature heat capacity C p on YM2B2C and LuNi2B2C in 
the vortex state have shown that C p clearly indicates the 
presence of the v^ff-term [Q . In the conventional s-wave 
superconductors, C p should increase linearly with H be- 
cause all the QPs are trapped within the vortex core with 
the radius of the coherence length £, and hence the QP 
density of states (DOS) N(H) is proportional to the num- 
ber of vortices; N(H) oc N F ^H where N F is the DOS at 
the Fermi level in the core ||. Thus, the observed y/~H- 
dependence of C p is strikingly in contrast to the ordinary 
s-wave superconductors. The -\/i?-dependence of C p has 
been reported in some of the unconventional supercon- 
ductors with gap nodes in the QP energy spectrum such 
as high-T c YBa 2 Cu 3 7 @ and heavy fermion UPt 3 |] 
superconductors. Nonlinear ii-dcpcndcnce close to the 



H of C p is also observed in some s-wave clean super- 
conductors such as CeRu2 M and NbSe2 [||||. 

Although the unusual ii-dependent C p of borocarbide 
superconductors has been discussed in terms of several 
intriguing models, the issue is still far from being settled. 
For instance, several authors proposed the shrinking of 
the vortex core with H |@,§-{y], but the physical ori- 



gin behind this phenomenon is unclear. Another group 
ascribed it to the field-induced gap nodes || , but this sce- 
nario is beyond the applicability of the original argument 
jjjj. Moreover, the extended QP states outside the vor- 
tex core owing to a presupposed <i-wave symmetry have 
been invoked [T^ |. There is, however, no corroborative 
evidence for d-wave pairing. Therefore, the experimen- 
tal clarification of the QP structure in the vortex state 
of YNi2B2C is very much needed. Further, it is crucial 
for understanding the vortex lattice structure. If a sub- 
stantial portion of the QPs extend in specific directions 
well outside the core, they should play an important role 
in determining the superconducting properties, including 
the vortex lattice structure, magnetization M(H), upper 
critical field H C 2, etc. Jlij . Nevertheless, these properties 
have been discussed in terms of an ordinary s-wave su- 
perconductor with an anisotropic Fermi surface with the 
use of the nonlocal London theory, without taking into 
account the effect of the extended QPs. Jljfl . 

We stress that we can extract detailed information 
about the QP spectrum only when we measure both the 
heat capacity and the surface impedance Z s for the same 
sample. In superconductors with gap nodes, C p is essen- 
tially determined by the QPs excited in the node direc- 
tions |fl6|| . Energy dissipation in the flux flow state, on 
the other hand, occurs mainly in the "normal regions" 
created by the vortices. Thus, the dissipative response is 
dominated by Andreev bound states localized within vor- 
tex cores; those states have momenta whose directions lie 
away from the nodes ]l7[ |. Therefore, the measurement of 
Z s is complementary to the heat capacity. In this Letter, 
we report measurements for both Z s and C p in the vor- 
tex state of YM2B2C. We have also examined the effect 
of introducing columnar defects (CD) on C p ; since the 
diameter of CD is comparable to £ ~ 70 A of YN12B2C 
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and the inside of CD is semiconducting p8| , they should 
strongly influence the electronic structure of the vortices. 
As we shall discuss below, these results provide strong ev- 



idence that the V-H-term in C p originates mainly from 
the Doppler shift of the delocalizcd QP spectrum due 
to the supcrfluid electrons surrounding the vortex cores 
p6[ . This suggests that a strongly anisotropic s-wave 
state most likely applies to YNi 2 B 2 C. 

Single crystals of YNi 2 B 2 C with T c = 13.4 K were 
grown by the floating zone method. In all measure- 
ments, dc magnetic fields H were applied parallel to the 
c axis. The microwave surface impedance Z s = R s + iX s , 
where R s and X s are the surface resistance and the reac- 
tance, respectively, were measured by the cavity pertur- 
bation technique. We used a cylindrical Cu cavity with 
Q ~ 3 x 10 4 operated at 28.5 GHz in the TE n mode. 
The sample was placed at the center of the cavity which 
is the antinode of the microwave magnetic field H ac that 
is parallel to the c-axis of the sample (H ac \\ H || c). In 
this configuration, the vortex lines along the c axis re- 
spond to oscillatory driving currents within the ab plane 
induced by H ac . We measured C p by the thermal re- 
laxation method. Several single crystals nearly one mil- 
limeter in size and ~ 100/im in thickness along the c axis 
were irradiated at GANIL (Caen, France) with a 6.0-GeV 
Pb-ion beam aligned parallel to the c-axis. The samples 
were irradiated to the fluence of 5 x 10 10 and 1 x 10 11 
ions/cm 2 , corresponding to a dose-equivalent flux density 
of -B$=l T and 2 T, respectively. The resulting dam- 
age consisted of a random array of amorphous columns, 
about 60-70A in diameter and continuous throughout 
the thickness of the crystal, which we confirmed by the 
transmission electron microscope image. The CD act as 
strong pinning centers of the vortices; at B < a 
substantial portion of the vortices are trapped inside the 
CD. Particularly below ~ 0.2B,/,, almost all vortices are 
expected to be trapped Jl^]. The CD do not, however, 
act as strong scattering centers of the electrons, unlike 
magnetic impurities or point defects. In fact, after the 
irradiation, M(H) showed a broad peak at B ~ B$ due 
to the trapping of the vortices by the CD (inset of Fig.l), 
but the increase in the resistivity was very small ( less 
than 10%) and the T c did not change. 

Figure 1 plots C p /T as a function of T 2 for the pris- 
tine YNi 2 B 2 C. We will first make a remark on the In- 
dependence of C p . Obviously, the normal state in Fig. 
1 (Cp/T vs T 2 ) that occurs as the superconductivity 
is suppressed by magnetic fields deviates from a linear 
dependence. Though this deviation may be attributed 
to the low-energy optical phonons, it is not clear where 
they originate from p0| . Moreover, recent inelastic neu- 
tron scattering experiments demonstrated an anomalous 
phonon branch which strongly couples to the electronic 
excitation in the superconducting state ^lj. Therefore, 
the quantitative analysis of the T-dependence of the elec- 
tronic contribution obtained by subtracting the phonon 
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FIG. 1. Cp/T versus T 2 for the pristine YNi 2 B 2 C. All mea- 
surements have been done in the field cooling conditions. In- 
set: M — H curves for pristine and irradiated =2 T) crys- 
tals at 5 K. While in the pristine crystal M(H) is reversible 
in a wide Jf-range, showing very weak pinning centers, M(H) 
show a broad peak at B = B<t, in the irradiated crystal. 

terms from the total heat capacity is precarious. 

Figure 2(a) depicts the -ff-dependence of the heat ca- 
pacity for the pristine YNi 2 B 2 C at a low temperature 
where the phonon contribution is negligible. In order to 
discuss the heat capacity induced by a magnetic field, 
we plotted AC P {H)/T = (C p (H)-C p (0))/T in Fig.2(a). 
A significant deviation from a linear ff-dependence is 
clearly observed. The inset of Fig. 2(a) plots AC P (H) /T 
as a function of \JlT . At a low field, AC P (H)/T increases 
with an upward curvature with respect to \f~Tl. However, 
at yH > 0.25 T 1 / 2 (H > 70 mT), AC P /T obviously in- 
creases linearly to This \/i7-dependence persists 
close to H C 2- Since the lower critical field H c i deter- 
mined by the M-H curve is of the order of 30 mT, the 
\Zff-dependent AC P /T is observed in almost the whole 
regime of the vortex state. Figure 2(b) shows the same 
plot for the irradiated YNi 2 B 2 C with B$=2 T. A sig- 
nificant deviation from the linear iJ-dependence persists 
even after the irradiation. Surprisingly, AC P /T is little 
affected by the introduction of the CD (see the inset of 
Fig. 2(b)). We obtained qualitatively similar results for 
YNi 2 B 2 C with Bg, — 1T. Before discussing these results, 
we will discuss Z s in the pristine YNi 2 B 2 C. 

The inset of Fig. 3 shows the T-dependence of Z s in the 
Meissner phase. Both R s and X s decrease rapidly with 
decreasing T below T c . Figure 3 shows Z s as a function 
of \H at 1.5 K. Let us quickly recall the behavior of Z s 
in type-II superconductors. In the Meissner phase, the 
microwave response is purely reactive and R s ~ and 
X s = [IowXl, where /io is the vacuum permeability, u> is 
the microwave frequency and Xl is the London penctra- 
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FIG. 2. (a) AC p (//)/T = {C p (//)-CV(0)}/T for the pris- 
tine YNi2B2C at low temperatures. Inset: AC P /T as a func- 
tion of y 7 //. (b) The same plot for the irradiated crystal with 
B* = 2 T @. Inset: The //-dependence of AC P /T before 
and after the irradiation. 

tion length. On the other hand, the response is dissipa- 
tive in the normal state and R s = X s = pqluS/2, where 
S = y // 2p n / hquj is the skin depth and p n is the normal 
state resistivity. We determined the absolute values of 
R s and X s from the comparison with the dc resistivity 
and assuming that R s = X s in the normal state. Using 
this procedure, we obtained Xl — 500 A. In the vortex 
state, Z s is determined by the vortex dynamics. Since the 
present microwave frequency is three orders of magnitude 
higher than the vortex pinning frequency of YN12B2C 
p2| , a completely free flux flow state is realized. In the 
flux flow state, two characteristic length scales, namely 
Xl and the flux flow skin depth Sf ~ y/2pf/Jiotj, ap- 
pear in accordance with the microwave field penetration. 
At a low field, A^ greatly exceeds Sf (Xl S> Sf). In 
this regime, R s and X s are given as R s ~ P//Xl and 
X s ~ hqujXl- On the other hand, at high fields where Sf 
greatly exceeds Xl {Sf ~> Xl), the viscous loss becomes 
dominant and the response is similar to the normal state 
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FIG. 3. Inset: T-dependence of Z a — R s + iX s in the 
Meissner phase. Main panel: R 3 and X s as a function of 
(/////c2) 1 ^ 2 at 1.5 K. The solid and dashed lines represent the 
result of the theoretical calculations of R„ and X 3 by Eq.(l), 
respectively, assuming two different //-dependences of N(H); 
N(H) cx H (case I) and N(H) cx VH (case II). For the detail, 
see text. 

(R s ~ X s ) except that 5 is replaced by Sf. In YN12B2C 
this crossover occurs at a very low field (poH ~ 20 mT). 
According to Coffey and Clem p3] ], Z s in the flux flow 
state is expressed as 



Z s = i^qujXl 



1 - ^/2)S f 2 /X\ 



•2iAi/<^ 



1/2 



(1) 



where S qp — y / 2p qp / pquj with the QP resistivity p qp 
is the normal-fluid skin depth. We analyzed the data 
by means of Eq.(l). We considered two different in- 
dependences of N(H), namely N(H) cx H (case I) and 
N(H) cx >/H (case II). Case I corresponds to the con- 
ventional Bardeen-Stephen relation, pf = (H / H r 2)p n 
and X\(H) = A|(0)/(1 - H/H c2 ), while case II corre- 
sponds to the vortex core shrinking, pf ~ \/H/ H C 2Pn 
and X L (H) = A|(0)/(1_- yjH/H c2 ). In both cases we 
used A L (0) = 500A |||. The solid and dashed lines 
represent the theoretical calculations of R s and X s by 
Eq.(l), respectively. Comparing the two cases, case I ob- 
viously describes the data much better than case II. Small 
deviations of the fits from the data in case I may be due 
to the fact that pf is not strictly linear in H, as suggested 
by the time-dependent Ginzburg-Landau theory |25f . 

The results of Z s and AC P for the pristine and irra- 
diated YNI2B2C offer important clues for understanding 
the QP structure. Figure 3 shows that the linear in- 
dependent Bardeen-Stephen relation yields a more con- 
sistent fit to the data than another model does. Since 
the QPs localized in the core mainly contribute to the 
flux flow dissipation, this fact shows that the number 
of the QPs trapped within each core is independent of 
in. Therefore, the scenario of the core shrinking with H 
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as an origin of nonlinear 7d-dependent C p proposed by 
several groups p|j9|-[Tl|] is completely excluded. In addi- 
tion and more importantly, the fact that the existence 
of the CD with a comparable radius with £ little affects 
the AC P implies that the QPs within the core radius £ 
are not important for the total heat capacity in the pris- 
tine YN12B2C. On the basis of these results, we are led to 
conclude that in YM2B2C the extended QP states around 
the vortex core play an important role in determining the 
superconducting properties, similar to d-wave supercon- 
ductors. 

Moreover, the effect of the CD on the AC P provides 
another important piece of information about the micro- 
scopic origin of v^ff-dependence. It has been pointed out 
that in the presence of line nodes there are two sources 
for the V^ff-dependent N(H), namely, the contributions 
from localized and delocalized fermions fl6[| , but which 
of the two dominates has thus far been left unquestioned. 
Here we shall attempt to address this issue. One arises 
from the localized QPs in the node directions. Since the 
extended QP wavefunction in the node directions is cut 
off by its adjacent vortices, the area per a single vortex 
is proportional to the intervortex distance R oc 1/ . 
Then N(H) is given as N(H) cx N F £RH cx The 
other arises from the delocalized QPs. In the presence of 
the supercurrent flow with velocity v s , the energy spec- 
trum of the delocalized QPs is shifted by the Doppler 
effect as E(p) — > E(p) + v s ■ p. In superconductors with 
line nodes such as d-wave symmetry where DOS has a 
linear energy dependence in the bulk (N(E) cx E), the 
Doppler-shifted QPs give rise to the finite DOS at the 
Fermi level. Then N(H) is obtained by integrating the 
vortex lattice cell N(H) cx R~ 2 v s (r) ■ p r dr cx \f~H, 
where |f s (r)| = h/2mr is the velocity of the circulat- 
ing supercurrents. The very weak influence of the CD 
on AC P elucidates the fact that the Doppler shift of the 
delocalized QPs is mainly responsible for the \/H depen- 
dent C p , because the contribution from the localized QPs 
in the node directions should diminish when the vortex 
lattice structure is destroyed. 

At low fields (H < 70 mT « H cr ), the AC P /T vs. V~H 
curve shows an upward curvature. There are two ori- 
gins for the deviation from the \f~H- behavior. The first 
one is the lower critical field H c \ ~ 30 mT, which is in 
the same order of H cr . The second one is the crossover 
from high field to low field scaling. Generally, the heat 
capacity of superconductors with line nodes is propor- 
tional to \IH only at high fields where the number of the 
Doppler-shifted QPs exceeds that of thermally excited 
QPs. The heat capacity at low fields should be described 
by the zero-field scaling C p /{^ n T) ~ fc^T/A, where 7„ 
and A are the Sommerfeld coefficient in the normal state 
and superconducting energy gap, respectively [p6| . The 
crossover field from high field scaling to low field scaling 
roughly estimated from the relation y/H/H C 2 ~ k-&T/ A 



with n H c2 = 5.5 T, T c = 13.4 K and T — 1.5 K is 
~24 mT, which is also in the same order of H cr . It 
should be noted that this behavior is also reported in 
YBa 2 Cu307_,5 with d-wave symmetry p7[ . 

The present experiments strongly suggest that the in- 
fluence of the extended QPs, which has been ignored by 
many authors, should be taken into account when dis- 
cussing the vortex lattice structure, H C 2 and M(H). In 
YM2B2C, de Haas- van Alphen oscillations are observed 
as low as H C 2/5 [ p8[ . This unusual phenomena may 
be related to the extended QPs. It is tempting to re- 
late the observed relevance of the delocalized QPs to a 
three-dimensional d-wave superconductivity |fL3f . How- 
ever, recent doping studies on YNi 2 B 2 C show that the 
superconductivity survives even in the dirty limit 0,^9) . 
This robustness of the superconductivity against impu- 
rities makes a d-wave state unlikely |30|. Therefore, we 
believe that a strongly anisotropic s-wave state in which 
N(E) cx E is similar to a d-wave state is most likely for 
YNi 2 B 2 C. 

We finally comment on the tunneling measurement by 
Ekino et al. which reports a full BCS gap in YN12B2C 
p3| . We believe that their break junction method selects 
a specific g-dircction in the gap function; in fact the same 
group reported a full BCS gap even in d-wave high-T c 
cuprates fl34| . 

In summary, we have studied the heat capacity and 
surface impedance of pristine and irradiated YM2B2C at 
low temperatures. We provided strong evidence that the 
Doppler shift of the delocalized QP spectrum around the 
vortex core plays an important role in determining the 
superconducting properties of YM2B2C. 

We thank N. Chikumoto, J. Clem, T. Hanaguri, 
R.P. Huebener, K. Maki, H. Takagi, A. Tanaka and 
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